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We predict the general feasibility and demonstrate the specific design of the THz laser operating 
between Landau levels in graphene placed on a polar substrate in a magnetic field of order 1 T. 
Steady state operation under a continuous wave optical pumping is possible due to an interplay 
between Auger and surface-phonon mediated relaxation of carriers. The scheme is scalable to other 
materials with massless Dirac fermions, for example surface states in 3D topological insulators such 
as Bi 2 Se 3 or BioTe 3 . 
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INTRODUCTION 


Free nonrelativistic electrons in a magnetic field behave as a system of harmonic oscillators, with selection rules 
allowing only the transitions between neighboring states with equal probabilities. Therefore they cannot be used 
as an active medium for lasers and masers. One way to get around this limitation is to accelerate electrons to high 
enough speeds that the relativistic effects become important. This leads to an anharmonicity in the electron spectrum 
and possibility of the maser action by accelerated electron beams, which has been so impressively implemented in 
vacuum electronic devices such as gyrotrons (Tj. Free carriers in semiconductors seem to offer a similar opportunity 
as the electron dispersion can show significant nonparabolicity above the bottom of the conduction band. Moreover, 
semiconductors offer a flexibility to grow heterostructures with different cyclotron transition energies which could be 
used for carrier injection into a given LL; see the proposal for a LL laser in the quantum Hall regime [2]. In practice, 
however, an ultrafast energy and momentum relaxation in semiconductors would quickly destroy population inversion 
between the Landau levels (LLs). As a result, there are no viable ’’solid state gyrotrons”, although Landau level 
quantization does help with reducing scattering rate and improving performance of quantum cascade lasers 0. that 
operate through population inversion between quantum well subbands. 

Graphene seems to be an ideal material for the realization of LL lasers. Low energy excitations near the Dirac 
points in graphene have a linear conical spectrum which is obviously extremely nonparabolic. In a transverse magnetic 
field the 2D conical spectrum splits into a series of non-equidistant LLs with energies scaling as a square root of the 
magnetic field and the principal quantum number. It was suggested in [3] that the optical pumping to an arbitrary 
excited state n > 1 will lead to electrons cascading down the LLs preferentially emitting photons, which would 
potentially lead to the EM field amplification on any of these downward transitions. Unfortunately, the proposal [3] 
assumed that the radiative transitions are the fastest ones in graphene. It did not include most important nonradiative 
relaxation channels and did not attempt to calculate actual LL populations. In particular, it turns out that the Auger 
relaxation is a very powerful relaxation mechanism for Dirac electrons in a magnetized graphene that proceeds much 
faster than radiative transitions and washes out any population inversion over the time scale of few ps; see below and 
also recent theoretical calculations of the Auger relaxation rate [Sj and experimental measurements in [B]. A recently 
proposed, more sophisticated pumping scheme j7| takes into account Auger relaxation processes and still leads to only 
a transient population inversion existing over a ps timescale. 

Here we propose what we believe is a viable inter-LL laser scheme for graphene that takes into account all relevant 
relaxation processes and in fact utilizes them to reach a steady-state population inversion, vital for any viable laser. 
Our scheme is transferable to thin (A A z > 5 nrn) films of 3D topological insulators such as BUSea where the 
Landau levels associated with massless metallic surface states BM should demonstrate similar coupling to the EM 
field despite different chirality (11)1 . Not that our scheme provides the population inversion in a steady state , i.e. under 
a continuous-wave pumping, in contrast to previous proposals, with or without the magnetic field, that could provide 
only a transient gain during a picosecond time interval mmm- 

We solve kinetic and density matrix equations coupled with Maxwell’s equations to calculate populations, gain and 
laser threshold conditions as a function of the optical pumping power. The calculation details are in the sections below. 
Here we present a general idea of the laser scheme. It is illustrated in Fig. 1. It shows one specific implementation of 
the scheme with an optical pumping originated from level n = — 2 to obtain maximum population inversion between 
levels -1 and -2. However, the scheme can me implemented for any pair of LLs (— n, —n — 1) as long as level —n stays 
deep enough below the Fermi level. The lasing wavelength can be from sub-THz to the mid-infrared range, depending 
on the value of n, the magnetic held, and the substrate used. 

The general idea is as follows. In the equilibrium (without pumping) labels n = — 1 and —2 are fully occupied 
to degeneracy surface density N s = g s g v /2irl c 2 , where spin and valley degeneracy factors are g s = 2 and g v = 2 
for graphene, and l c = \Jcfi/eB is the magnetic length. The Fermi level is placed at the Dirac point in the figure, 
assuming intrinsic graphene. However, this can be changed, as long as level n = — 1 is fully occupied in equilibrium. 
An optical pumping resonant to the transition —2 —> 1 moves part of the carriers up from level n = —2 creating a 
population inversion between a fully occupied level n = —1 and level —2. In order for this population inversion to 
exist in a steady state, i.e. under a continuous-wave optical pumping, the relaxation of carriers back to lower laser 
state —2 should be slower than the relaxation rate to the upper laser state —1. Unfortunately, the Auger mechanism 
does not satisfy this criterion. Our simulations show that there is no population inversion in the steady state, no 
matter how strong the optical pumping power is. This is because an increase in the depopulation rate of level —2 
by an optical pumping is compensated by an increase in the Auger scattering rate to level 2, primarily through the 
scattering of electrons from states in levels 1 and —1 to states in levels 2 and —2, respectively. In order to overcome 
this obstacle, the magnetic field needs to be tuned in order to bring the transitions 1 —>■ 0 and 0 — > — 1 (of the same 
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FIG. 1. The scheme to obtain population inversion between the electron states below the Fermi level by using a continuous-wave 
optical pumping. 


energy) in resonance with an LO phonon energy. This will greatly increase the rate of electron relaxation from excited 
states to the upper laser state n = — 1 through LO phonon emission, whereas the transitions to state —2 will be out 
of resonance and not affected much. 

The LO phonon energy in graphene is close to 200 meV, which would require a magnetic field of almost 30 T to 
bring the transition frequency wio close to ujlo- In order to reduce the required magnetic field one can utilize the 
scattering on bulk, surface, or interface optical phonons of the substrate, and choose the substrate with a lower optical 
phonon energy, for example a polar semiconductor such as GaAs or InGaAs |T2]. For definiteness, below we assume 
the substrate to be GaAs, which leads to the surface optical (SO) phonon energy of hcogo = 36 meV [IB] . This is 
equal to uiqi = ui c = '/2vf/l c in a magnetic field of 1 T. The laser transition wavelength would be then around 82 
//m, i.e. around 3 THz, which is the range where there is a shortage of laser sources. We will also assume the optical 
pumping between levels —2 and 1, although the pumping resonant to the transition from —2 to 3 would be equally 
efficient and lead to a similar value for the gain. Moreover, the transition frequency for the latter transition in a 
magnetic field of 1 T would correspond to a CO 2 laser wavelength around 10 p m, which could be more convenient 
than the 14-/mr wavelength corresponding to the transition —2 —» 1. Of course all energies can be changed as needed 
by choosing different substrates or different LLs for the lower laser state, for example n = —3 instead of —2. 


ELECTRON STATES AND OPTICAL TRANSITIONS BETWEEN THE LANDAU LEVELS IN 

GRAPHENE 

For completeness, we give a brief summary of the electron states and optical transitions between the LLs in graphene, 
since this information is extensively used below. They have been calculated many times before and observed both in 
monolayer and multilayer samples El- 

Neglecting intervalley scattering, we will only need electron states in one of the two equivalent K, K' valleys, for 
example the K valley. Without a magnetic field, the low-energy Hamiltonian in the vicinity of the K Dirac point is 
given by m 


H = vf<? • p = vf 


0 Px - iPy 
Px + iPy 0 


(1) 


where Vp = 10 8 cm/s. In the presence of a transverse magnetic field or any EM field described by the vector-potential 
A, we replace p with n = p + eA/c. For a magnetic field in the +z direction, we can write A = (0, Bx , 0) in the 
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Landau gauge, then the eigenfunctions are expressed as m 

F«(f) = -^e ik y<f> n (k,x), (2) 

with 


®n{k,x) = C n 


/sgn(ra)iH 1 <f>\ n \_ 1 (x + %k)\ 

V i^4>\n\{x + lc k ) ) 


( 3 ) 


where C n = 1 when n = 0, and C n = l/y/2 when n/0; sgn(x) = 1, 0, -1 for x > 0, x = 0, x < 0 respectively. The 
function (j)\ n \ (x) has the same form as the eigenfunction in the massive electron case: 


0|n| (*^) 


1 

— , - exp 

^/2\ n \\n\\y/Trl c 


1 

2 




( 4 ) 


where H\ n \{x) is the Hermite polynomial. The corresponding eigenenergy is e n = sgn(n)hoj c y/\n\, with oj c = %/2 vf/i c . 

In this manuscript, an electron state will be labeled by |n, fc, s, £), where s = {f, 4-} denotes spin, £ = { K , K'} 
denotes valley; fc, s, £ are degenerate quantum numbers, and the total degeneracy density of a Landau level n is 2 / 7 W 2 . 
The interaction Hamiltonian for an optical field with an in-plane polarization can be written as 


HZ = v F - c a-A°v 


( 5 ) 


where A op is the vector potential of the optical field, which is related to the electric field by E op = (—l/c)dA op /dt. 
If we define two circular polarization vectors, l ® = (x + iy)/y/ 2 and Iq = (x — iy)/y/2, the vector potential of a single 
frequency optical held can be written as 


A° p = i (Aj e + A e l e ) e~ iwt + c.c. . 


( 6 ) 


Plugging this expression into the Schrodinger equation and using the rotating wave approximation, we get the same 
selection rules as in m-- Transitions between rq and 712 [n-i > n{) are coupled by photons with 1 0 polarization if |ri 2 | 
= \n\\ + 1 , and with Iq polarization if |ri 2 1 = |tii| - 1 . 

By expressing A op through E op in HZi we can § e ^ the magnitude of the dipole moment for a resonant transition 
between Landau levels ni and U 2 - 


|Mnin2| — y/^C ni C n2 CV F /iO . 


( 7 ) 


The two dimensional linear optical susceptibility near the resonance to the transition between n\ and ri 2 ( n\ < 712 ) is 


Xnirt 2 — 


\l^r, 


\fn 2 - fn J 


nl 2 c huj - (e„ 2 - e ni ) + ih/T 2 
where T 2 is the dephasing time. The optical transition rate between n\ and n 2 (n 1 < n 2 ) is 


F°p = 

mn 2 2 


V J n 1 n 2 E op 


1 /T 2 


(1/T 2 ) 2 + ((e„ 2 - e ni )/h - w) 2 


( 8 ) 


(9) 


LASER THRESHOLD CONDITION 

To determine the threshold condition for a LL graphene laser we consider the simplest geometry resembling a 
quantum-well vertical cavity laser, in which an active layer consisting of one or several graphene monolayers on a 
polar substrate is located between the two mirrors of given reflection factors see Fig. 2. We will assume that 
there are two media with dielectric constants and k 2 from both sides of the active layer. We will also assume for 
simplicity that the thickness of an active layer is much smaller than the wavelength of the THz laser field. For a field 
of amplitude Ei incident on the graphene layer, the amplitudes of reflected and transmitted waves E r and E t can be 
related using the Maxwell’s equations with proper boundary conditions as 
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FIG. 2. A design of graphene laser. 
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( 10 ) 


X 


To make equations even simpler, we will take K\ = n 2 =k. It is straightforward to include more complex cavity 
structures if needed for a particular design. 

The fields also need to satisfy the boundary conditions at the mirrors: 

E u e~ ikL ' = n{E lr + E 2t )e ikLl , 

E 2i e~ ikL2 = r 2 (Ei t + E 2r )e ikL2 . (11) 

From the boundary conditions Eqs. (10) and the condition to have stable nonzero optical fields inside the cavity 


is 


27 xiuj 


-X = 


T\v 2 — e~ 


( 12 ) 


J‘^T 2 -\- V\ 6 21fcZ.2 _J_ i*2e _j_ g 2ik{L\-\-L2) 

To get the threshold, we assume that the optical fields are in resonance with respective transitions and the lengths L\ 
and L 2 are adjusted so that every term in the denominator has the same sign. Then the minimum required imaginary 


part of the susceptibility in the active layer can be found from the real part of Eq. (12): 

kir 2 | 


=^Mx] = yt i i 
V KC i + |n| 


l - 


(13) 


|?"2| + l r l r 2| 

We will discuss the feasibility of reaching the lasing threshold below, after calculating the rates of scattering processes, 
the non-equilibrium populations of the LLs, and the resulting graphene susceptibility at the laser transition in the 
presence of an optical pumping. 

The carriers excited by an optical pumping relax through a variety of scattering processes. The steady state 
populations are determined by a balance between relaxation and the continuous wave pumping. In the next two 
sections we give a detailed description of most important processes that determine the redistribution of populations 
and the resulting steady-state gain. 


AUGER PROCESSES 

A strong magnetic field suppresses scattering processes due to energy quantization and reduction in the phase space 
available for scattered carriers. However, Auger processes remain very efficient: due to the symmetry between electron 
and hole LLs there is always resonance for scattering of carriers from (0,0) LLs into (1,-1) states and for all other 
combinations allowed by the energy conservation: (1,-1) o (2,-2), (0,0) O (2,-2), (—1,1) o (2,-2), etc.; see 
Fig. 3. Recently the Auger relaxation rates were measured to be in a few ps range in pump-probe experiments [6], 
which agrees with our simulations. Below we outline the general derivation of the Auger scattering rate and then 
apply it to our problem. 
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FIG. 3. Examples of the Auger scattering processes between n = 0, ±1, ±2 LLs. . 


General formulas 


Auger processes are mediated by the Coulomb interaction between carriers. The general Coulomb interaction 
Hamiltonian for electrons can be written as JT5] 




1 

2 




(14) 


where 


v afhS = wmmivcouiin - f 2 )| 7 (i))i(j(2)). (is) 

In order to simplify this expression and include the effect of screening, we expand Vcw(ei — ri) in Fourier series 

Vcoui (n - f 2 ) = X] V ? e<f ( * - ( 16 ) 

Q 

where V$ = 2ne 2 / KoAq for a 2-dimensional case. Using this expression, we get 


Vafrt = ^^(a(l)|e^|7(l))</3(2)|e-^|5(2)) . 
<? 


(17) 


To include screening, we replace V r j with V s (q,ui) = V$/e(q,ui), where the dielectric function e(q,ui) in the random 
phase approximation is given by the Lindhard formula 

e(g,w) = 1 - F f n°(5,a;) , (18) 

and the polarizability n°(<f, w) is written as 


n°(£ W ) = £ — 


f a ~ Si 3 


a/3 


ep + hu ; + iS 




(19) 
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with the form factor F a p(q) = (a\e lq ' r \fi). The value of w is determined by hcu = E 1 — E a in V a p 7 s '1E|- The rate of 
the Auger scattering from state |a, b) to state |c, d) is calculated from the Fermi’s golden rule; it is symmetric with 
respect to the initial and final states: 


Tabbed = \{cd\V c \ab)\ 2 5{E C + E d - E a - E b ) , 


( 20 ) 


where the matrix element is 


(cd| ^ (y^edab Vdcab H - Vdcba Vcdba) • 


( 21 ) 


So, there are essentially four terms because electrons are indistinguishable particles. The state | a,b) in Eq. (21) only 
means that both |a) and | b) are occupied, instead of specifying that electron 1 is in |a), and electron 2 is in | b). One 
can find mistakes in the literature with some of the terms missing. 


Auger scattering between Landau levels in graphene 


For graphene in a transverse magnetic field, an electron state can be written as |a) = |?r Q , k a , s a ,t; a ), with notations 


explained in Sec. II. We will address the screening effect first. The form factor in Eq. (19) can be evaluated to be 


F a p(q) = (ale^l^) = <J aai8fl ^ ai f^, #ifca _^e-^ , c^G nanp ( gi ,,« I! ) . 
So the polarizability becomes 


( 22 ) 


n / \ 2 a 4 v —a 

n (fo W )= W2 E 


fn a ~ f, 


np 


7r// c —' e„ 

np 


- e np +hco + iS 


\Gn a np idlyi 9x) | 


(23) 


One can check that \G nanfl (q y , q x )\ only depends on the magnitude of q. so II 0 (g, ui) and V s (q,u>) are only functions of 
q = \q\. We calculated the form factor numerically and checked that it agreed with the analytical expression in [T5] , 
Then, using Eq. © and Eq. ( |29| ) , the Coulomb matrix element is 

Vabcd = ^Sa,,S c ^Sb,Sd^a,ic^b,id 

x ^ v s&u)5 qyika -k c e- iq * l2 c k *G nana (q y ,q x )5- qy ,k b -k d e iq:Deokd G nbrid (-qy,-q x ) 


^ S a.,s c fiSb ,s d<5£ a ,£c^£b lid^ka + kb ,fc c +fcd 


q y =k a — k c 


(24) 


For a fixed k a , this matrix element decays quickly when k c — k a is large, since G nanc (q y , q x ) oc exp(—(gZ c ) 2 /4) when q 
is large m- If k c is bounded, then k b is bounded too, otherwise the term e lq - r ^d k <- k d) WO uld oscillate too fast with 
q x , which essentially makes the summation over q x to vanish. 

The Auger scattering rate between two pairs of Landau levels ( n a ,n b ) and ( n c ,rid ) is 


1 


n a nb<->n c rid 


2 A/nP c 


E E E E T ab<r+cd ( x - if n a = n b or n c = n d 


^ci^ci^c ^di^di^d 


(25) 


where the factor 1/2 in the parenthesis is because of the double counting the initial or final states. In Eq. (251, one of 
the summations can be dropped immediately, since the result from the other three summations will be independent 
of the forth set of quantum numbers. This summation will give exactly the degeneracy 2A/ttI x , so it will cancel with 
the pre-factor. One summation of k can also be eliminated due to the conservation of momentum k a + k b = k c + k d . 
As the energy is fully quantized, we will replace the S function in Eq. (20) with a Lorentzian of line width which can 
be attributed to impurity scattering m- 










PHONON SCATTERING 


General formulas 


The interaction Hamiltonian between phonons and electrons can generally be written as 

H£t = E F ^) c \ + r^+ b -*> ’ ( 26 ) 

k, q 

where c and c 1 are annihilation and creation operators for electrons, b and E are annihilation and creation operators 
for phonons, and F(q) is defined below. Using Fermi’s golden rule, the scattering rate from an initial electronic state 
| tpi) to a final state | <pf) is 

r E = ^ EK + l)\F(q)\ 2 \M fi (q)\ 2 S(e f + Ef - €i ) 

Q 

o?r 

+f Y,n q \F{ q )\*\M fi {q)\ 2 ^f-EZ h - e 4 ) , (27) 

q 

where the first term is for the phonon emission, the second term is for the phonon absorption, and the matrix element 
is given by 

M fi {q) = (</?/lE4 +Fk\Fi) 

k 

= J drip* f (r)e l? " r ipi(r) (28) 

Using the wave functions in Eq. ([ 2 ]), the matrix element Mf i(q) for |z) = |nj, fcj, s, £) and |/) = |n/, kf, s, £) is calculated 
to be 

S qy , kf - ki e~ iqxl - ki G nfni {q y ,q x ) , (29) 

where we have defined 

G ni n 2 (q y ,qx) = J dx^l 1 (q y ,x)e tqxX ^ na (0,x) . (30) 

The averaged scattering rate from an initial Landau level rii to a a final Landau level n/ is 



E Y'Ph 
k t 

Y EK + l)\F{q)\ 2 \G nfn My, qx)\ 2 5(e f + E* - e,) 

Q 

^Y 2 nq\F(q)\ 2 \Gn f nMv,q x )\ 2 S(e f - Ef - Cj ) . 

Q 


(31) 


We include the dynamic screening effect by carriers in graphene in the phonon scattering processes. This can be 
done by replacing F(q) with F s (q,ui) = F(q)/e(q,u>), where the dielectric function e(q,u>) is given in Eq. (18) and 
Eq. (23). We will only consider the case of low enough temperatures, when phonon absorption is unimportant and 
only phonon emission processes contribute to the scattering rate. At room temperature this is still a reasonable 
approximation; it can be easily dropped if a greater accuracy is needed. 


LA phonon scattering 


For longitudinal acoustic (LA) phonon scattering, the expression for F(q) is [201 


Fla(q) 



D^/q , 


(32) 
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where p = 7.6 x 10 -8 g/cm 2 is the area mass density of graphene, v s 
in the 10 - 50 eV range. Also, the energy of a LA phonon is E^ A = 
we get the scattering rate by LA phonons: 


= 2 x 10 6 cm/s is the sound velocity, and D is 
hv s q. Plugging these expressions into Eq. (311, 


■pZ/.A 


D 2 2 r2ir 

47r ^° 2 J d8\G nf ni (q 0 sin9, q 0 cos6)\ 2 , 


(33) 


where go = (e ni — e nf )/hv s . The coefhcient is of the order of 10 14 s 1 for B 


1 T. However, the integrand in Eq. (33) 
= a/2 Vf/v s = 50-/2. So, 


is roughly of the order of exp[— (go^c) 2 /2], which is extremely small, since qol c ~ u c l c /v s = V2 vf/v s = 50/2. So, the 
LA phonon scattering does not contribute significantly to electronic transitions between Landau levels due to a large 
ratio Vf/v s 1 . 


Surface optical phonon scattering 


Since we want to use phonon scattering to our advantage in order to facilitate electron relaxation to the upper laser 
state, we consider graphene on a polar substrate or sandwiched between two substrates. In this case, the electrons 
in graphene can couple to the surface or interface modes of optical phonons pa Eij, which we will call the surface 
optical (SO) phonons for brevity. If the two substrates on both sides of the graphene layer are the same, the SO 
phonon energy is equal to the longitudinal optical (LO) phonon energy of the substrate [2TJ . If there is vacuum on one 
side, the SO phonon energy is slightly shifted from the LO phonon energy m- We will assume the former case for 
definiteness, but note that it would be straightforward to calculate interface optical phonon modes for an arbitrarily 
complex structure. If we assume that the graphene layer does not affect the SO phonon modes, then the expression 
of F(q) can be written as 


Fso{q) 


2 'Ke 2 hu>so 

f 1 

1 \] 

A 

x-sub 
\ ^00 

^sub I 

/ J 



(34) 


where A is the area of graphene, Kq (k^°) is the low (high) frequency dielectric constant of the substrate, and huso 
is the energy of the surface optical phonon. Since it has a flat dispersion, we replace the S functions with a Lorentzian 
Cj(E) = 7 / 7 t(E 2 +7 2 ), where 7 is the broadening of Landau levels, which can be again attributed to disorder. Using 
again Eq. (31), we find the SO phonon scattering rate to be 


r SO _ 1 2 , 
f USO 


^sub ^sub 

^00 ^0 


•£7 €rtf 


dq 


0 ‘(9-^110(5,0,)) 


/ d6\G nfni (qsin 9, qcos0)\ 2 , 

Jo 


(35) 


where the screening effect is included, and ui = (e ni — e nf )/h. 


LANDAU LEVEL POPULATIONS UNDER OPTICAL PUMPING 


After the expressions for the optical transition rates and the scattering rates due to SO phonon emission and 
Auger processes have been found, we can write the density matrix equations with adiabatically eliminated optical 
polarizations to arrive at the set of rate equations for the filling factors of the Landau levels: 


d 

dt 



op d 

SO 7 

a 

+ dt tn * 

+ jt Ua 


Auger 


(36) 


where 


d 

dt 


op 


A 

dt 


SO 


Y, T Zn b (fn a ~ fn b ) , 


Y, T n°^nJn a {i - f nb ) + ^rr^/^l - /„J , 


n b 


n b 


(37) 


(38) 
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and 


~f 

dt ln 


Auger 


— ^ v ^ a ,fi 


*E E r. 

n c n d >n 


rib 

n a nb<->ri c nd 


(~fnjn b { 1 - /nj(l - fn d ) + fnjn d ( 1 - /nj(l - fn J) 


(39) 


where t 


n a ,rib 


= 2 if rib = n a , and 1 otherwise. 


Using these rate equations, we can simulate the dynamics of the graphene system for an arbitrary optical excitation. 
Note that the system is highly nonlinear, firstly because of the state filling and secondly, because the matrix elements 
depend on the dynamic screening, which depends in turn on the instantaneous distribution of electrons in Landau 
levels. Therefore, time dependent simulations are time consuming. Here we present the steady state results for the 
continuous-wave optical pumping. 


RESULTS AND DISCUSSION 


For a GaAs substrate, the SO phonon energy is 36meV, which requires the magnetic field to be around IT. In 
the simulations, broadenings of all transitions are set to be 5meV, and T 2 is 0.1 ps. Also, we consider intrinsic 
(undoped) graphene as an example, so without pumping the n = 0 LL is half-filled, all LLs below are fully filled, 
and all LLs above are empty. We define the gain between n = — 1 and —2 as the left-hand side of Eq. (131: g. lt . 2 
= —(27ruj/-^Kc)Im[x(u>-i,-2 )]■ To minimize the absorption of the THz field by the polar substrate, we would like to 
reduce its thickness to a few fim to be much smaller than the wavelength of the THz field but at the same time, thick 
enough to be considered bulk for SO phonon scattering. The gain is maximized when there is air outside the active 
layer so that k = 1. 

The dependence of the steady state filling factors and gain per graphene monolayer on the pump intensity are 
shown in Fig. [4] and Fig.[5j 
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FIG. 4. Dependence of the steady state filling factors on the pumping intensity. 


As can be seen from the figures, one can achieve a significant steady-state population inversion between states 
n = — 1 and —2 and the gain value of about 0.05 per monolayer of graphene. One can scale the gain up by stacking 


many graphene monolayers. For comparison, the right-hand side of Eq. (131 which describes mirror losses is equal 


to 0.025 when the reflectivities rq = = 0.95, which is easily achievable. The closest allowed transition at the Iq 
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FIG. 5. Dependence of the gain between n = -1 and -2 LLs per graphene monolayer on the pumping intensity. 


polarization is from n = —3 to n = —2 LLs. It is detuned from the laser transition frequency by about 4 meV in a 
magnetic field of 1 T. Therefore its contribution to losses is lower than the gain. Since the electron motion is quantized, 
there are no other losses in graphene associated with free carriers. The undoped GaAs is a popular material for the 
nonlinear THz generation and its THz losses are rather low, especially since the polar substrate can be thinned down 
to a few fim. Therefore, one can operate in the desirable regime where the losses are dominated by mirror losses. 

For surface states in 3D topological insulators such as Bi2Se3 or Bi2Te3, the Fermi velocity has a similar value but 
there is no spin or valley degeneracy. Therefore, a similar laser scheme with a thin film of Bi 2 Se3 (i.e. two surfaces) 
placed on a polar substrate will give about two times smaller gain. Additional free-carrier THz losses may exist in 
this case due to unintentional doping of the bulk Bi2Se3. 

One can also see in Fig.[4]that the population inversion exists also between states n = 2 and 1, albeit at a two times 
lower level. This seems unexpected, given that the optical pumping brings carriers only to state 1. However, a closer 
look at the rate equations shows that the population inversion between levels 1 and 2 is a consequence of a strongly 
non-equilibrium carrier distribution below the Fermi level created by the optical pumping, namely the population 
inversion between states —2 and —1. Indeed, when f_ 1 > /_ 2 , the Auger scattering rate from states (1, —1) to states 
(2, —2) is greater than the scattering rate in the opposite direction. This creates the population inversion / 2 > /i and 
the gain for the l® polarization, which is about two times smaller than the Iq gain. 

In conclusion, we show the feasibility of the Landau level THz laser in a magnetized graphene. Despite ultrafast 
Auger relaxation, steady-state operation of the laser under continuous wave optical pumping is possible by utilizing 
surface or interface phonon relaxation. The scheme is scalable to thin films of 3D topological insulators such as Bi 2 Se3 
or Bi 2 Te 3 . 
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